Optical conductivity of one-dimensional doped Hubbard-Mott insulator 
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We study the optical response of a strongly correlated electron system near the metal-insulator 
transition using a mapping to the sine-Gordon model. With semiclassical quantization, the spectral 
weight is distributed between a Drude peak and absorption lines due to breathers. We calculate the 
Drude weight, the optical gap, and the lineshape of breather absorption. 



Optical absorption is one of the basic tools in stud- 
ies of the spectral properties of strongly correlated elec- 
tron systemsaBn. An important class of such systems 
is comprised of one-dimensional conductors such as lad- 
der compounds, organic metals and stripes in high T c 
cuprates. The electron-electroa interaction in these low- 
dimensional systems is knownB to result in spin-charge 
separation and formation of a Luttinger liquid or the 
Mott-Hubbard insulator (MHI). The MHI occurs at com- 
mensurate filling fractions and is characterised by an en- 
ergy gap in the charge sector. In this case, optical absorp- 
tion exhibits a threshhold at twice the gap energy. As the 
Fermi level shifts above the gap, the system conducts, 
with profound changes to the optical absorption spec- 
trum. The optical conductivity summarizes this behav- 
ior: <t{lo) — D 2t:S(lu) + cr rcg (w), with the Drude weight 
D — in the insulating phase, and a regular part a con- 
taining contributions from pair excitations and bound 
states, or excitonstl 

The simplest model where the MHI phase occurs is 
the Hubbard model with an on-site repulsion. In the 
commensurate phase the spectrum of the optical absorp- 
tion consists of- |[; ap | broad feature lying in the particle- 
hole continuurroErO. In the more realistic extended Hub- 
bard model, the excitation spectrum generally consists 
of both thej-Harticle-hole continuum and bound states 



(cxcitons)liM Numerical studies of the extended Hub- 
bard model with nearest- and next-nearest neighbor 
repulsionalj show the emergence of a sharp exciton peak, 
which becomes a dominant feature when the repulsion is 
sufficiently strong. 

In addition to the numerical analysis of the problem a 
mapping to the continuous sine-Gordon (SG) model can 
be considered!! The exciton peaks are seen in the sine- 
Gordon approach as absorption on quantum "breather" 
modes of the model. 

Much is known about the commensurate phase be- 
cause of the integrable structure of the corresponding 
SG theory and the knowledge of its form-factors due to 
the Smirnov bootstrap (seeEj and refs therein). The in- 
commensurate phase is much less studied with analytic 
results existing for the high frequency expansion of the 
optical conductivity in the absence of breathersa. 

In this paper we report our results for the optical con- 
ductivity in the incommensurate phase in the presence 
of a large number of breathers, which corresponds to 
long-range interactions in the Hubbard chainOJ. It is 



known that in.-this limit the SG model can be analyzed 
semiclassicallytM Using the semiclassical approach we 
calculate the optical absorption of a doped MHI as a 
function of frequency and the carrier density. Our ana- 
lytic expression for the Drude weight (|2l| ) is found in an 
excellent agreement with a thermodynamic Bethe ansatz 
calculation. At high doping the optical weight is con- 
tained in the Drude peak and in a series of well resolved 
peaks at finite frequencies, corresponding to excitation of 
breathers. We calculate the positions of the peaks jis|), 
and their shape ( p3| ) . At low doping peaks merge in a nar- 
row continuous band and we provide an analytic expres- 
sion for the optical absorption in this case. Our results 
for one spatial dimension agree with experimental and 
theoretical evidences for the peaks in optical absorption 
at small energies near the metal-insulator traflsijtic 



quasi ID compounds and in higher dimensionso ljEllijtLj. 

The optical response, o£ a doped MHI is described by 
the sine-Gordon modelfl'Eil 

£ = - 2 {(dt<j>) 2 ~ (d x ^) 2 } + cos/3^ + h^-dA (1) 



where the charge and current densities read 

P = -a.*, 3 = -^. 



(2) 



In Eq. (jlj) h is the chemical potential controlling the aver- 
age density (p(x)) = p, and it is assumed that in the insu- 
lating phase p = 0. The cosine term in Eq.([l]) is related to 
Umklapp. The sine-Gordon couplinsu/3 is determined by 
the interactions of the lattice modeltl For free fermions 
f3 = \/8n. For infinite U Hubbard model /3 = V47T, 
which corresponds to the refermionization point, where 
the optical absorption can be found exactlyEjj. Below, 
we will investigate the case of (3 <C 1 when the number 
of breather modes is ~ 8w/f3 2 . he mass m is a charac- 
teristic energy scale for a breather excitation while the 
particle- hole continuum starts from energies ~ 8m//3 2 . 

By virtue of @) the Kubo formula for the optical con- 
ductivity is written as 



cf(lj) = — Im / dte 
\2nJ Jo 



where 



G{t)= / dx([ir(x,t),<f>(0M)- 



(3) 



(4) 



2 



and 7r = dt4> is the canonically conjugate momentum for 
(f>, such that [<j)(x), n(y)] = i5(x — y). 

It follows from (0) that <r(b>) satisfies a sum rule 



duj cr(b)) = i 



An 



dx{[w(x,Q),<j>(a,Q)]) = 2t, (5) 



with t = /3 2 /87r. In the Tomonaga-Luttingcr limit, m = 
in (|l|) , the sum rule @ is saturated by the Drude weight 
D = t/tt. 

In the quasiclassical regime the case p = is very sim- 
ple. Eq. (|l|) reduces to the Klein-Gordon model, with 
the ground state <p = 0. The optical conductivity is cal- 
culated immediately 



(tkg(w) = r (8(b) — m) + 8(b) + m)) , 



(6) 



For a finite p the quasiclassical ground state found from 
the stationary SG equation under constraint (d x <f>) = 
2irp/f3 is 



(j) (x) = 2(3- 1 &m(2Kpx) 



(7) 



where am(i, k) is the Jacobi amplitude functiorJlS with 
the elliptic index k. This index is found from the equation 



m/p = 2kK(k) 



(8) 



with the complete elliptic integral K. Henceforth we will 
drop the index k in the argument of the elliptic func- 
tions and use the conventions k\ = y/l — k 2 , K{k) = K, 
K(k\) = K' , E(k) = E. The standard semiclassical ap- 
proach requires the second variation of the action around 
the classical solution and an analysis of the piectrum of 
the resulting Gaussian action. It was, showrJla that this 
scheme encounters infrared divergences in higher orders 
of perturbation theory in /3. An infrared stable formula- 
tion of the theory is achieved by a non-linear change of 
variables 



(f> = (f)Q (x + k 1 (3rj) . 



(9) 



where k = 2-np. The new field variable t](x, t) relays the 
role of the continuous "collective coordinate" iB, which 
roughly corresponds to fluctuations in kinks' positions. 
The linearized Lagrangian for the field rj reads 



C = \w\x){(d t r 1 f-(d x r,f}. 

where the weight function, w is defined as 

. , 2K , 
w(x) = an2Kpx. 



(10) 



(11) 



The spectrum of the Lagrangian ( |10| ) was analyzed in de- 
tail inta. The eigenmodes of (nG) are the Bloch functions 



T)q(x,t) 



iqx—ibJ q t 



Xg(x), 



(12) 



where q is the Bloch wave vector and x is the modulating 
Bloch function with period 1/p. 
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FIG. 1: Semiclassical result for the spectrum in the reduced 
Brillouin zone scheme for p/m = 0.14. The first three allowed 
transitions are labeled wi ,W2 and u>3 . 



These solutions are conveniently parametrized by a 
complex parameter a as follows. The Bloch wave vec- 
tor and the eigenfrequency of the mode are given by 

q = -i2pKZ(a) + np, lu = 2pKdn(a), (13) 

where Z(a) is the Jacobi Zeta function. The modulating 
Bloch function is 

Xq (x) = C a e^x (irpx - ^) ^\npx). (14) 

Here the coefficient C a is chosen to fulfill the normaliza- 
tion condition 



dxw 2 (x)r]a(x)r]p(x) = 8 a p. 



(15) 



The quantized field rj(x,t) is then represented in the 
oscillator basis 



/2uT r 



h.c. 



(16) 



with operators b, satisfying usual commutation rela- 
tions [b a ,bl,] = 5 aa >. 

In the limit of small f3 the change of variables (|J) can 
be linearized 4> = 00 + W7 li n = W V anc l horn the nor- 
mal mode expansion (|l^) and the Kubo formula (^) one 
obtains for the optical conductivity 

a(b)) =Tj2\F{a)\ 2 {5{u;-Lu a ) + S(LU + u) a )) 1 (17) 

a 

where the form factor 

F(a) = -)= ( dxr] a (x)w(x). (18) 
vL Jo 



The selection rules for the form factors in (17) should 
allow only transitions between the states of zero wave 
vector in the reduced Brillouin zone scheme (see Fig. [l]). 
This condition combined with ( |l3| ) implies that the al- 
lowed indices a and the corresponding frequencies of the 
absorption lines satisfy 



2iKZ(a n ) — ■n{2n — 1), b) n — (m/fc)dna rf 



(19) 



3 



where n = 1, 2, 3, . 
factors we find 



and a n S [0,iK'). For the form 



\F(a n )\ 2 = 



tt 2 csch 2 [i^+iK'Z(a n )] 



2K 2 



dn 2 av, 



E_ 

K 



(20) 



The lowest frequency io\ can be associated with the opti- 
cal gap. From Eq. ( |l9|) (see also Fig. [j]) one finds that uj\ 
is always greater than m and is an increasing function of 
p. One has ui\ ~ m for p <C m and wi ~ 2irp for p~^> m. 
The spacings between the peaks also increase with p, 
from zero at p — > to 2irp for p ^> m. The amplitudes 
|-F(a„)| 2 are rapidly decreasing functions of n. In the high 
density limit one finds |-F(a n )| 2 ~ n~ 2 (4Tip/m)~ An . 

The Drude weight comes from the zero mode 770 and is 
found to be 



7T 1 1 4EK 



(21) 



Combined with definition (|J) formula ( pl| ) gives the qua- 
siclassical expression for the Drude weight as a function 
of the density of kinks. The following limiting cases allow 
for further simplifications: 

1. High density limit. This limit corresponds to k — > 
and K,E — > tt/2. For the Drude weight we get D — ► t/7t, 
which is the result for the Tomonaga-Luttinger limit. 

2. Low; density limit. In this limit both A: and E are close 
to unity. The Drude weight is then given by 



D 



2AL 



(22) 



This result has a simple physical interpretation. In the 
limit of low density the kinks form a gas of weakly inter- 
acting non-relativistic particles. In an external electric 
field e the motion of a particle is described by Newton's 
law —iM s cuv = e and the electric current is given by 
j = —pe/(i<jjM s ). Analytic continuation of this expres- 
sion gives Eq. (|22|). 

To check the accuracy of these results we calculate 
the exact Drude weight from the thermodynamic Bethe 
ansatzEl A comparison made in Fig. [2] shows an excel- 
lent agreement between the quasiclassics and the exact 
result. 

Further analysis needs a different treatment of the 
cases of the high and the low density. 

In the limit p/m — > one has K — > 00 so that both the 
weight of 5— functions in ([l?]) and the spacings between 
them vanish. (In this limit one has K' — > 7r/2, E — * 1, 
Z(a) — > tanha, dna — > secha and a n — —i(2n+l)Trp/m, 
uj 2 = in 2 + {(2n + l)7rp] 2 . ) As a result, the sequence of 
5— functions merges to the continuous function 



mp 



d(uj 2 — m 2 ) 



>\Vu 



m 2 sinh 2 -f—^JuJ 1 

2m v 



(23) 



Note the non-integrable singularity ct(uj) ~ (to — m) 3 / 2 
at the threshold, which seemingly violates the sum rule 
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FIG. 2: Drude weight as a function of soliton density in the 
small t limit. The solid line corresponds to the semiclassical 
result, Eqs. (§), @. The exact values of the Drude weight 
calculated from the Bethe ansatz for r = 0.05 are shown by 
circles. 



(|5j). This signals that the low density limit requires spe- 
cial treatment. A careful analysis of the fluctuations in 
this limitcil leads to the formula 



o-(w) 



dxe 



1£-x coth x 



(24) 

which satisfies the sum rule and coincides with ( J23| ) for 
LU — 771 > p 2 /m. 

Next we turn to the case of large densities. The series 
of infinitely narrow peaks in Eq. ( |r7| ) stems from the clas- 
sical picture of the broken translational symmetry of the 
ground state Eq. (Q) which is a periodic lattice of kinks. 
The true ground state, however, should be translation- 
ally invariant, for arbitrarily small (3. The symmetry is 
restored by zero-point fluctuations of the respective po- 
sitions of kinks in the lattice. Below we show that such 
fluctuations lead to non-perturbative corrections to Eq. 
( |l"7| ) and to the algebraic broadening of the conductance 
peaks. 

In order to analyze finite (3 corrections to Eq. ([l7]) it is 
convenient to write the following identity for the correla- 
tion function M) 



(25) 



G(t)= j dx—([ S mP<j>(x,t),<j>(0)]), 



which follows from the equations of motion of the sine- 
Gordon model. The optical sum rule is taken care of by 
imposing the initial condition G(0) = — i. 

With a help of (|) and observing that the ground state 
solution (j)Q satisfies the stationary SG equation we write 



G(t) = / dx{[$ \ x + 



Prj{x,t) 



/5??(0) 



/x 



(26) 



where the average over the position of the lattice of kinks 
X is spelled out explicitly. The function O is a sum of 
a linear and a periodic function of its argument and can 



4 



be represented as a Fourier series 



o \x + 



f3r]\ kx 



+ ^ — e il( - KX+Pv) (27) 



k J (3 



with Ci = sech(TrlK'/K). Using @ we write Eq. @ 



as 



G(t)= K 2 L- 1 f dxdy^Cfe^-^Q^x - y,t), (28) 



where 



Mz-y,*) = /r 2 ([e 



ilf3ri(x,t) e -iZ/3»j(2/)l 



(29) 



In deriving Eq. (pq), off-diagonal terms of the form 
([exp(im/3?7(x, i)), exp(— in/3ri(0))]) for n ^ m were 
dropped, since they vanish in the thermodynamic limit. 
The average over the position of the kink lattice X was 
replaced by the equivalent integration over y. 

The limit f3 — > in (p8|) is non-trivial because of the 
long-distance singularities characteristic of this problem. 
First take this limit naively by expanding the exponen- 
tials in Eq. ( |29[ ) in powers of small (3 and keeping the 
leading term. Then 



Qi{x-y,t)nl 2 {[ V {x, t),ri(y)}). 



(30) 



Using this approximation in formula ( p8| ) immediately 
leads to the peaked structure of the optical response (17). 

Next, consider the average commutator in Eq. ( |29| ) 
using the Gaussian approximation 



, e i/3l V {x,t) e -if3l v (0)\ 



(31) 



Eq. ([30]) is obtained from the first two terms of the 
Taylor expansion of the exponential in the right hand 
side of (|3l|). This expansion is accurate while the argu- 
ment of the exponential is small, which is not true at 
large distances, where p 2 (rj(x, t)rj(0)) « /Clog \x 2 — v 2 t 2 \, 



where K. is the Luttinger parameter,t2l. The breakdown 
of the approximation Eq. ( |30|) at large distances leads to 
a broadening of the conductance peaks. 

A rather lengthy but straightforward analysis, based 
on Eq. (|3l]) shows then that the real part of optical 
conductivity near the n-th peak is given by 



a(uj n + il) cx R- c 



dt 



l 2 C 2 e~ mt 



[p 2 {v 2 n - v*)t 2 + 



. (32) 



where v n is the group velocity of the bosonic excitations 
near frequency u n , v s is the sound velocity of bosonic 
excitations. We can see, that in the limit JC = the 
delta-function form of the conductivity peak is restored. 
At any finite value of the Luttinger parameter K, the con- 
ductivity will fall off algebraically away from u> n . 

a{uj n + ft) cx l^l 2 ^ 1 (33) 
The width of the peak, on the other hand, is determined 
by the difference of the sound velocity v s and the group 
velocity v n and is given by 



AO = pJv 2 - vl 



(34) 



In the high density limit one finds Afl m 2 / p which im- 
plies, that peaks of optical conductivity are well resolved 
despite the quantum smearing of the lattice structure. 

In summary, the optical conductivity of a one- 
dimensional conductor close to the Hubbard-Mott transi- 
tion in the quasiclassical limit shows a smooth crossover 
between the Luttinger liquid and the Mott-Hubbard in- 
sulator. With increasing density the band of breather 
absorption splits into a series of well resolved peaks. We 
calculated the positions of the peaks (|l9|), their shape 
( p3| ) and width (|34|). Our analytic result for the Drude 
weight (21) shows an excellent agreement with the ther- 
modynamic Bethe ansatz. 

We thank A. Rosch and P. Horsch for helpful discus- 
sions. 
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